Bound whispering gallery modes in circular arrays of dielectric spherical particles 
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Low-dimensional ordered arrays of optical elements can possess bound modes having an extremely 
high quality factor. Typically, these arrays consist of metal elements which have significantly high 
light absorption thus restricting performance. In this paper we address the following question: 
can bound modes be formed in dielectric systems where the absorption of light is negligible? Our 
investigation of circular arrays of spherical particles shows that (1) high quality modes in an array 
of 10 or more particles can be attained at least for a refractive index rir > 2, so optical materials 
like Ti02 or GaAs can be used; (2) the most bound modes have nearly transverse polarization 
perpendicular to the circular plane; (3) in a particularly interesting case of Ti02 particles (rutile 
phase, rir = 2.7), the quality factor of the most bound mode increases almost by an order of 
magnitude with the addition of 10 extra particles, while for particles made of GaAs the quality 
factor increases by almost two orders of magnitude with the addition of ten extra particles. We 
hope that this preliminary study will stimulate experimental investigations of bound modes in low- 
dimensional arrays of dielectric particles. 

PACS numbers: 61.43.Fs, 42.25.Fx, 42.55. f, 71.55.J 



1. Assemblies of microsize particles are useful in a 
variety of optical applications because of their resonant 
interaction with a visible light. Although most of the 
present research is focused on three dimensional (3 — d) 
systems (e. g. photonic crystals), low-dimensional (1 — d 
and 2 — d) structures are attracting increasing attention 
because they can be more easily constructed and main- 
tained than their 3 — d counterparts. Despite its low 
dimension, 1 — d structures have a remarkable array ge- 
ometry effect on absorption, photoluminescence and Ra- 
man scattering of lighti^ii^ They can even serve as a 
1 — d nano- waveguide^ 



The strong dependence of light scattering and absorp- 
tion on 1 — 0? array geometry suggests that the quasi- 
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FIG. 1: Circular array of particles. Arrows indicate polar- 
izations of whispering gallery modes, tl along the z-axis; 12 
along the radii and I tangential to the circle, both in the x — y 
plane. 



bound optical modes exist in extended I — d arrays. The 
resonant interaction of light with those modes is respon- 
sible for the optical properties of arrays. Each mode 
can be characterized by the complex eigenfrequency z = 
Lu — ij/2, with the real part uj representing the mode fre- 
quency and the imaginary part 7 representing the mode 
decay rate.^ Modes of interest must have a large enough 
quality factor Q — uj/"f so that they are bound to the 
whole 1-d structure. 

Indeed, 1 — d arrays of N identical particles can possess 
bound modes with a quality factor approaching infinity 
for N oo,SiL2i2iia so -f{N = 00) = 0. The case Q ^ 00 
takes place in a linear chain of identical particles sepa- 
rated by a small distance a that is less than half of the 
resonant wavelength A 

a < A/2 = c/(2w). (1) 

If Eq. is satisfied, then there exists at least one quasi- 
bound mode having a maximum wavevector qmax = 7r/a 
exceeding the wavevector of the resonant photon k = 
2tt/X. When the chain is infinite the decay of such mode 
is forbidden by the momentum conservation law requiring 
k > q. In the finite system of N particles, the decay rate 
is finite, but it tends to zero when N ^ 00 in accordance 
with the power law ^{N) oc N^^ (see Ref(&). 

Quasi-bound modes of the finite system possess a very 
narrow frequency resonance that can be used in a vari- 
ety of photonics applications including waveguides)^ an- 
tennas and detectors liSiS These applications involve the 
emission or absorption of light in a very narrow frequency 
range, and the capability of functioning in the single pho- 
ton regime.'-^ This narrow resonance leads to a well re- 
solved far field emission pattern that can be used as a 
guiding signal for aircraft. "'^^ In addition, systems pos- 
sessing modes with a high quality factor can be used in 

lasers because these modes have a low pump threshold 
for lasing.13,14^15,16 
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2. One can expect that the highest quahty factor can 
be attained for whispering gaUery modes in an array of 
particles arranged in a circle and equidistant from each 
other (Fig. ^ . A circle has no sharp ends in contrast with 
a particle chain where the lifetime of a mode is limited 
to the travel time between the ends of the chain. ^ In- 
deed, it has been demonstrated that whispering gallery 
modes in a circular array of cylinder-shaped antennas 
perpendicular to the circle plane possess a quality fac- 
tor that grows exponentially as the number of anten- 
nas increases iiSiiLiS) Similarly, as described in Ref»ii, 
an exceptionally high quality factor of circularly shaped 
structural fluctuations was suggested to account for las- 
ing in random media. Note that these modes can also 
produce a large Raman surface enhancement.!^ 

We believe that the exponentially small decay rate of 
whispering gallery modes is a common property of ap- 
propriately constructed circular arrays. This expectation 
can be justified as follows. Consider an array of N iden- 
tical particles forming a circle of radius R and separated 
by distance a (see Fig. QJ. The symmetry of the prob- 
lem suggests that rotation by 2tt/N radians will change 
eigenmode amplitudes by a factor of exp{iqa), where the 
wavevector q = m/R {m = 0, — 1, 1, — 2, 2, ..iV/2), enu- 
merates A'" eigenmodes possessing different quasi-angular 
momenta m. We will consider only even numbers of N for 
the sake of simplicity. The mode can decay by the emis- 
sion of a photon. The angular momentum of the emitted 
photon takes on the discrete set of values rn, m -\- 
m — N, ... due to rotational symmetry. We are primarily 
interested in the most efficient regime, when the angular 
momentum is equal to its minimum value m. The rate 
of emission of a photon possessing angular momentum m 
is defined by the squared overlap integral of the bound 
guiding mode and the photon wavefunction. This inte- 
gral can be estimated using the wavefunction of a photon 
within the circular array (Fig. ^ positioned as character- 
ized by the cylinder coordinate p — R. This wavefunction 
is given by the Bessel function Jm{kR), where k — to/c 
is the wavevector of an emitted photon, and lu is its fre- 
quency. Thus one can estimate the whispering gallery 
mode decay rate dependence on the number of particles 
N in the circular array as 

j{k)^\JgR{kR)\^,N = 27TR/a. (2) 

This result reproduces the estimate of RefiiS for interact- 
ing antennas in the very similar case of a mode possessing 
a certain quasi-angular momentum. If the argument of 
the Bessel function is smaller than its index {k < q in ac- 
cordance with Eq. Q), then using the limit N ^ oo one 
can approximate Eq. |(2J) by the exponential functioni^ 
(see Fig. 

"/ (X exp{—KaN/Tr), (3) 
K — q{cosh^^{l/x) — \/l — x^), X = k/q. 
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FIG. 2; Quality factor of g = vr/a modes versus the number 
of particles for tl, and / modes in Ti02 (rutile phase, rir = 
2.7) and ZnO and t\ mode in GaAs {ur = 3.5). The fit 
of the quality factor by an inverse Bessel function Eq. Q 
Q = 0.2598/ I JkaN/(2TT) P is shown for the longitudinal mode 
in Ti02 (the product qa has been defined by the solution of 
Eq. mill at A = 110. The mode t2 is not shown because it is 
hardly distinguishable with the mode tl at A S> 1. 



Thus our estimate shows that the radiative decay rate 
of the whispering gallery mode can be made exponen- 
tially small. Exponentially narrow resonances (Eq. ^) 
can contribute to the optical properties of the system 
only when the width of the frequency resonance (Eq. Q) 
exceeds the rate of light absorption by particles. Metal 
particles or antennas always have a significant absorp- 
tion due to conducting electrons. This is not the case 
for dielectric particles. They have negligible absorption 
when the photon energy is less than the width of the for- 
bidden band. In addition, the optical (Mie) resonance 
frequency of a dielectric particle is inversely proportional 
to its size and can be easily changed to the desired value, 
while the resonance of metal particles is bound to their 
plasma frequency at least until the particle size is less 
than the wavelength. Therefore we find it important to 
test whether narrow resonances (Eq. can be attained 
in the array of dielectric particles possessing a realistic 
refractive index n^. 

3. Bound modes have to satisfy Eq. Formally, 
this goal can be attained by reducing the interparticle 
distance a. However, the interparticle distance cannot 
be made smaller than the spherical particle diameter d, 
so we cannot satisfy the condition k < q a,t any refrac- 
tive index. Since the resonant wavelength increases with 
increasing refractive index Ur approximately as A oc n^, 
bound modes will definitely exist at sufficiently large n^. 
Standard optical materials have refractive indices of or- 
der of 1 (see Table IJ). It is not clear whether arrays 
of particles made of those materials can possess bound 
modes. Below, we employ the multiple sphere Mie scat- 
tering formalism.^^ (see also2i*2^) generalized to the eigen- 
mode problem, and show that bound longitudinal modes 
exist at least for particles made of GaAs and Ti02 and 
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for any other material with the refractive index Ur > 1.9. 
The bound transverse modes should exist at arbitrary 
refractive index at sufficiently large number of particles 
as one can expect in analogy with the weakly guiding 
optical fibers 

Eigenmodes of the array can be found by solving corre- 
sponding Maxwell equations without incident field. One 
can express electric and magnetic fields at a given fre- 
quency z (wavevector k = z/c) by expanding them over 
spherical vector functions weighted with scattering am- 
plitudes a^„, for each sphere I {I — 1,2,3, ...N), 
where index n — 1, .. stands for the photon angular mo- 
mentum and index m = —n, —n + 1, ..n is the angular 
momentum projection onto the z-axis^ Eigenmode am- 
plitudes a and b are bound by the set of equations 

= ^ + (B<„ + (lC„, (4) 

where coefficients a[j, are the Mie scattering coeffi- 
cients of the sphered In particular the dipolar scat- 
tering coefficients oT, bi are defined as^i 

_ _ ■)jji{kd/2)ip[{nrkd/2) -nripi{nrkd/2)-)p[{kd/2) 
^ c^i{kd/2)i;[{nrkd/2) - nrMnrkd/2)c:[{kd/2) ' 

— _ nripiikd/2)ip[{nrkd/2) - ipi{nrkd/2)i>[{kd/2) 
^ ~ nr';i{kd/2)i(;[{nrkd/2) - il;i{nrkd/2)c;[{kd/2) ' ^ ' 
ipi{x) — sin[x)/x — cos{x),(,i{x) — e'"(— 1 — i/x). 

Matrices A and B define the interaction of multi- 
pole polarizations of different spheres. Their matrix el- 
ements have a general structure -M^„^j^ = exp{i{m — 
")'^ii)Sp-Dp(^ii)e2;p(i^r/j)/rf-, where nj, cj)ij, 9ij are 
the spherical coordinates of the center of the sphere j 
with respect to the center of the sphere I; k = z/c; and 
Dp{9) are the real coefficients as defined in RefiSS. The 
solution of Eq. Q exists for the discrete set of frequen- 
cies Za — (^a + ija, a — 1, 2, 3.., whcrc these frequencies 
are defined by eigenfrequencies of modes uja and their 
decay rates 7a. 

The rotational symmetry of the problem (Fig. ^1 per- 
mits us to seek solutions of Eq. I^J in the form 

(n\ — n p'2Tri{m+qa)l W _ i 2iTi(m+qa)l ij-s 

with wavevector q = 2TTp/{Na) {p = -1, 1, -2, 2. ..N/2). 
The dimension of Eq. Q is reduced by a factor of N, 
because the intersphere interactions A and B get re- 
placed with their Fourier transforms. At the next step 
one should restrict the maximum value of the angular 
momentum to n < rimax and solve the equations requir- 
ing the weak sensitivity of the solution to the increase 
of Umax- In this paper, we will study the simplest case 
n-max = 1 for coefficients bmn, while we set all coefficients 
o-mn equal to zero. This approach is similar to that of the 
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FIG. 3: Dependence of the quality factor of longitudinal mode 
on the refractive index near its threshold value n* ~ 1.9 esti- 
mated for = 100 of dielectric spheres within the array 

coupled dipolar approach. It is well justified for a large 
refractive index ^ 1 and for modes corresponding to 
dipolar Mie resonances possessing the lowest frequency^^ 
In fact, we are interested in these modes because they 
are most easily bound (cf. Eq. 

Optical materials like Ti02 or ZnO (see Table IJ) do 
not possess a high enough refractive index to use this 
dipolar approach. However, bound modes as found in 
the dipolar approach will remain bound as higher multi- 
poles are taken into account. In fact, higher multipoles 
add additional modes as generated by corresponding Mie 
resonances, and all are located above the dipolar reso- 
nance. Those resonances can shift the energy of dipo- 
lar resonance downwards only because of energy level re- 
pulsion. Such behavior has been seen, for instance, in 
Refj^, where the rigorous approach is compared with the 
discrete dipolar approach. The approximating method 
underestimates the size of the sphere at the lowest Mie 
resonance, which is equivalent to overestimating the reso- 
nant frequency. Therefore Eq. will be better satisfied 
in the more rigorous approach because of the increase of 
resonant wavelength due to the reduction of resonant fre- 
quency. Thus our method gives upper estimates for the 
decay rate of bound modes and for the minimum refrac- 
tive index n* necessary to form strongly bound modes. 



In the dipolar approach, one can separate whispering 
gallery modes into three branches tl, t2 and / depending 
on their polarization (Fig. The particle polarization 
for the transverse mode tl is perpendicular to the array 
on the X — y plane and is described by the amplitude 
floi (m = 0). The transverse mode t2 is characterized by 
particle polarizations close to the radial directions of the 
X — y plane. It is defined by amplitudes a±ii (a_ii w 
2aii-- ). The longitudinal mode I is characterized by the 
polarization that is approximately tangential to the circle 
and is defined by amplitudes a'^-n (a_ii w — 2aii). In the 
case of most bound mode with the angular momentum 
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q = 7r/a, the definitions of transverse and longitudinal 
modes are precise. 

The Fourier transform of interactions in the simplified 
Eq. (0J) using the parameter definitions of Ref.^» yields 
the dispersion equations for all three modes ijl = tl, t2 
and Z, given by 

1 

=-r^ +ef^{k,q) ^0,k = z/c, (8) 
ai[z) 

characterized by the dispersion laws 

etiik, q)^^ (Si(fc, q) + iE2(fc, q) - ^sik, q)) , (9) 



et2{k,q) - -{e+{k,q) + ^ e^{k, qf + V {k, qY (10) 

ei{k,q) = \{e+{k,q) - ^ e^{k, qf + V {k, q^ 

, , (t(/c. (7 + 27r/a) ± cr(fc, g — 27r/a) 
^±\k,q) = ^ ; 

a{k, q) = 0.75(Ei(A:, q) - i^2{k, q) + ^^{K ?)); 
V{k, q) = 0.75(Si(fc, q) + 3^E2(fc, q) ~ S^siK q)). (11) 

where Ep(g) = X^zIT^ ea;p(ifcro/ - iqal)/rli, tqi = 
2Rsin{7Tl/N) is the distance between the centers of 
spheres separated by Z — 1 spheres (see Fig. and 
q± = q zL 2-k/N . For N oo, the dispersions of modes 
tl and t2 become identical because the circular array of 
large radius locally approaches a linear chain geometry 
where two transverse polarizations are equivalent. The 
"photon" field term Tii(k,q) disappears in the Z— mode 
spectrum Eq. Hill) for iV ^ oo because photons have a 
transverse polarization.^'' 

4. One can resolve Eqs. lO, ^(TU\i . ((TT|l using the 
Newton-Raphson iteration algorithm for the equation 

fik)^^}-r + e^{k,q)^0 (12) 

with the iteration procedure defined as 

fc„+i fc„ - /(fc„)//'(fc„). (13) 

This procedure starts with some initial value k — kQ. At 
almost any initial value ko, the iteration series rapidly 
converges to some solution of Eq. (|12|l . However, Eq. 
(|12|l has generally infinitely many solutions and we are 
only interested in the solution possessing the largest qual- 
ity factor. In particular, if there exists a quasi-bound 
mode possessing a quality factor much greater than unity, 
then the procedure should be organized in such a manner 
that it converges to this mode. In the regime of interest 
of Eq. (Q, where the bound mode exists, we have orga- 
nized the iteration procedure choosing the starting point 
ko = 2/d. Using this choice, the algorithm converges to 
a frequency of the bound mode with the smallest decay 
rate in almost 100% of our studies for tl modes. The 
situation is more difficult for t2 and I modes, where, at 



certain number of particles the algorithm Eq. H13|l does 
not converge being trapped by the attractor formed by 
the selfrepeating sequence of two or three different values 
of g appearing periodically in Eq. H13|l . In these cases the 
convergence has been reached by changing initial value 
of q to 0.2/ d or using the modified Newton Raphson al- 
gorithm fc„+i = kn — 0.1/(fc„)//'(fc„). When the quality 
factor is high {Q > 20) this is sufficient to find the right 
solution. 

In our calculations, we have used the particle diame- 
ter d = 2 and took the minimum possible interparticle 
distance a — d, where bound modes can be formed most 
easily (Eq. |^). Although the accurate analysis of elec- 
tric and magnetic fields will be difficult in this case due 
to the singularity at the point where particles touch each 
other, the dipolar approach should still serve as a good 
upper estimate for mode frequencies and decay rates as 
discussed above. The generalization to different particle 
sizes is straghtforward because characteristic frequences 
and wavectors behave as g, fc, z oc l/d, while the quality 
factor is scale-invariant. The resonant frequencies for the 
particle size d = 200nm for a large number of particles 
A'^ = 10 having different polarizations and refractive in- 
dices are given in Tabled for the case of the most strong 
bound modes with q = tt/o. We have also included mode 
quality factors for A^ = 10 and their rates of exponential 
increase with the number of particles A^. Fig. H illus- 
trates the exponential dependence of the quality factor 
on the number of particles for different modes and mate- 
rials. 

Indeed, for large number of particles (A^ > 20), the 
smallest decay rate has always been found at the max- 
imum wavevector q — ■n/a corresponding to the quasi- 
angular momentum m — N/2 (cf&SiS and we further dis- 
cuss only this case). As transverse modes have the low- 
est frequencies uii{nr,N), these modes are the most eas- 
ily bound. Bound longitudinal modes are formed when 
uji{nr,q = vr/a) < ujthr = cn/a. According to our cal- 
culations, this happens at > n* ~ 1.9. The behav- 
ior of the quality factor of longitudinal mode near the 
threshold is shown in Fig. O For Ur < 1.9 the quality 
factor for a circular array of A^ = 100 dielectric particles 
depends very weakly on the refractive index, while the 
sharp, stepwise increase is clearly seen when the thresh- 
old is passed. Qualitatively similar behavior takes place 
at larger A^. Also for Ur < 1.9, the quality factor is almost 
independent of the number of particles (cf. Fig. ^ while 
above this point the dependence is quite strong. As for 
transverse modes they can be made bound at arbitrary 
refractive index^ However for rtj. < ?^*, a very large 
number of particles is required to attain a high quality 
factor (see Fig. El for ZnO). 



5. Thus we can conclude that strongly bound lon- 
gitudinal and transverse whispering gallery modes exist 
in circular arrays of Ti02 and GaAs particles because 
their refractive indices exceed the upper estimate of the 
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Material 


ZnO 


Ti02, anatase 


Ti02, rutile 


GaAs 


rir 


1.7 


2.55 


2.7 


3.5 


uJt (s-') 


2.2 ■ 10^^ 


1.8 ■ 10^^ 


1.6 • lO^'' 


1.2 • lO^'' 




2.6 ■ 10^^ 


1.9 ■ 10^^ 


1.8 • 10^^ 


1.4 ■ 10^^ 


Q, iV = 10, tl 


12.7 


86.77 


124 


828 


Q,N = 10, ; 


7.8 


40.3 


55 


303 


dln{Q)/dN, tl 


0.0157 


0.175 


0.21 


0.4 


dln{Q)/dN, I 





0.099 


0.131 


0.308 



TABLE I: Characterization of most bound modes in circular arrays of particles of 200 nm diameter. Mode parameters have 
been calculated for A'^ = 10 particles. 
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FIG. 4: Contour plot describing the longitudinal whispering 
gallery mode intensity for A = 10 GaAs particles {rir = 3.5) 
circular array. The intensity is expressed in arbitrarily units 
and contour lines separate domains for the intensity change by 
factor of ten (increase or decrease of the decimal logarithm of 
intensity by one. Particle sizes here and in Fig. |S]are chosen 
d = 200nm. 



threshold n* = 1.9 (see Fig. OJ. Quahty factors of most 
bound modes at the maximum wavevector q = n/a for 
ZnO, Ti02 and GaAs are shown in Fig. |2| They aU 
are exponentially sensitive to the number N of particles 
and one can relate this dependence to that of the inverse 
squared Bessel function Eq. ||2Jl as illustrated in Fig. [21 
for TiOa. 

To characterize the spatial structure of whispering 
gallery modes, we have created contour plots for the in- 
tensity distribution of bound longitudinal modes within 
the array x~y plane for GaAs particles. The most bound 
modes {q = 7r/a) are shown in Figs. 0I[31for = 10 and 
A^ = 20 particles respectively. In both cases the mode 
energy is highly concentrated within the array and the 
shift from the array domain by the distance comparable 
or less to the particle size leads to the reduction of the 



FIG. 5: Contour plot describing the longitudinal whispering 
gallery mode intensity for A = 20 GaAs particles (n^ = 3.5). 



mode intensity by orders of magnitude. The localization 
of optical energy is clearly stronger in the case A^ = 20. 
Indeed, relatively small displacement by 300nm from the 
domain of array in any direction leads to a reduction of 
intensity by the factor of 10^ (only factor of 10 in the 
case of A^ = 10). This behavior agrees with the sharp in- 
crease of the quality factor with the number of particles 
for GaAs (see Fig. 12). The strong reduction of intensity 
near the center of the circle is the consequence of the large 
angular momentum of the mode under consideration (cf. 
Eq. in the limit r ^ 0). 

As was discussed, the actual frequencies of low energy 
modes are smaller than our estimate, so the decay rates 
should decrease faster with A^ than predicted (see Eq. 
(|5|)). Even within our approach, the decay rate of the lon- 
gitudinal mode for = 3.5 (GaAs) decreases by almost 
two orders of magnitude with the addition of ten particles 
to the array. To attain a similar effect in a 3 — d photonic 
crystal, it can require the addition of 10'^ 1, 000 parti- 
cles, which is much more difficult. Therefore, we believe 
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that remarkable progress in photonics technology can be 
made by using 1 — d circular arrays of dielectric particles. 
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